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Abstract 

The boundstate problem in 2+1-dimensional large-A^ Yang-Mills theory is 
accurately solved using the light-front Hamiltonian of transverse lattice gauge 
theory. We conduct a thorough investigation of the space of couplings on 
coarse lattices, finding a single renormalised trajectory on which Poincare 
symmetries are enhanced in boundstate solutions. Augmented by existing 
data from finite- Euclidean lattice simulations, we obtain accurate estimates 
of the low- lying glueball spectrum at N = oo. 
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I. INTRODUCTION 



Yang- Mills theories are theoretically interesting in 2 + 1 dimensions because their prop- 
erties are very similar to the corresponding theory in 3 + 1 dimensions, yet they can be 
handled much more accurately; see Ref. for a review of properties and extensive refer- 
ences. They appear to exhibit linear confinement of heavy sources, a discrete spectrum of 
(glueball) boundstates, and a finite-temperature transition. Teper has recently performed a 
comprehensive analysis, using the standard tools of Euclidean SU{N) lattice gauge theory, 
of 2+1-dimensional Yang-Mills theories at = 2, 3, 4, and 5. Hamiltonian lattice calcula- 
tions have also recently been performed for finite and, though less comprehensive, 
the results are mainly consistent. With data at enough values of A^, one can contemplate an 
extrapolation to = oo. This is a limit of special interest for 'analytic' approaches to gauge 
theory, which often take advantage of large- A^ simplifications. In the absence of any other 
criteria for the errors involved, the only way to know how well 'analytic' approaches are do- 
ing, for example those of Refs. [H,^, is to compare with lattice data and their extrapolation 
to large A^. 

A related question is: how close is A^ = oo to small A^? This question can only be 
faithfully answered once there are accurate results in both limits. The 1/A^ expansion is 
typically an asymptotic one and, a priori, observables in the two limits need not be close in 
value. The existing finite-A^ data suggest strong suppression of corrections to the large- A^ 
limit [0,0, a conclusion that was speculated about much earlier on the basis of less 
reliable lattice data. If true, this fact deserves a deeper understanding. 

The main objective of this paper is to address these issues for 2 + 1 dimensional Yang- 



Mills theory with explicit calculations at A^ = oo. In Refs. (see also |jTj]) we used the 

large- A^ limit of 2 + 1-dimensional Yang-Mills as a test for developing the transverse lattice 
method of solving non-abelian gauge theories [0. Based on an improved understanding of 
the sources of error in that calculation, we perform here a calculation at the next level of 
approximation. We obtain a renormalised light-front Hamiltonian on the transverse lattice 
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for both the pure-glue and heavy-source sector. From this, we obtain the glueball spectrum 
and the heavy-source potential. Existing finite- data, combined with our explicit large- 
results, are used to determined the first few coefficients of the l/N"^ expansion of glueball 
masses in string tension units. 

In the next section we briefiy review the transverse lattice method, the details of which 



have been covered elsewhere [0,|r2,|l^]. Section 3 describes the numerical search for the 
renormalised Hamiltonian via tests of Poincare invariance. Our thorough investigation yields 
a single, well-defined candidate for the renormalised trajectory in coupling space. Results 
for the low-lying glueball eigenstates on this trajectory and the first few coefficients of the 
l/N"^ expansion of their masses are given. This improves upon current estimates of the 
large- limit, allowing us to accurately verify that corrections to it are highly suppressed. 
In the conclusions we discuss possible reasons for the success of analytic approaches, given 
their approximations. 

II. TRANSVERSE LATTICE IN 2 + 1 DIMENSIONS 

Adapted to 2 -|- 1 Yang-Mills theory, the Bardeen-Pearson transverse lattice gauge the- 
ory consists of continuum gauge potentials {^0,^2} and space-time co-ordinates 
together with gauge- covariant transverse link variables M{x^) running between sites at 
and -|- a on a transverse lattice of spacing a. We also use the light-front combinations 
x± = (x° ± x2)/V2, = (^0 ± A^)/V2, et cetera. DLCQ pf and Tamm-Dancoff cut- 
offs on the number of partons are used as intermediate regulators. These are extrapolated 
following the analysis of ref. [|ri[|. DLCQ means that we impose anti-periodic boundary 



conditions x~ ~ x~ + 27rfC/P"^, where is the total light-front momentum, and K is an 
integer cut-off. remains continuous and infinite, and is used as a canonical time vari- 
able to derive a light-front Hamiltonian P^. The most general action, from which P~ will 
be derived canonically, must allow all gauge invariant operators that respect the Poincare 
symmetries unviolated by the (gauge invariant) cut-offs. Since we will explicitly extrapolate 
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the DLCQ and Tamm-Dancoff cut-offs, only local dimension 2 operators with respect to 
{x~^,x~} co-ordinates will be included at the outset. In the discussion hereafter we assume 
this limit has be taken. 

Near the transverse continuum limit a ^ corresponding to SU{N) Yang-Mills, one 
expects M to take values in SU{N). However, away from this limit one can allow M to 
be a general NxN complex matrix, provided it still gauge transforms covariantly. One 
must then search this larger class of lattice theories for the renormalised trajectory that 
leads one to the continuum limit SU{N) Yang-Mills theory. Physical results are invariant 
along this trajectory and equal to the values in the full continuum limit. The trajectory 
may be found by renormalisation group transformations in the neighborhood of a fixed point 
(continuum limit). However, this is difficult for the present formulation. There are (roughly) 
two possibilities for the behaviour of M at a given point in the space of couplings constants: 
M is a massive degree of freedom (M = is the classical minimum); or, the 'radial' part of 
M condenses. We expect the latter to be the case near the a = limit of Yang-Mills, where 
the action should be minimized near values of M in SU{N) rather than M = 0. Dealing 
with the condensation of the radial part, or using unitary matrices for M from the outset 
0], is tricky in light-front quantisation. This is what makes an analysis near a = difficult. 
On the other hand, it is straightforward to perform canonical light-front quantisation about 
M = 0, when this is a stable minimum. If the renormalised trajectory passes into such a 
region, we can then study it. 

An alternative way to find the renormalised trajectory is to use symmetry |16|. Generally 



speaking, we can define a quantum field theory by symmetry — in our case gauge and 
Poincare invariance — and a particular continuum limit (there may be more than one). 
There is actually no reason why we cannot take a partial continuum limit, a limit in some 
space-time directions but not in others, since Poincare invariance should relate them. Thus, 
in Ref. |]1T| we proposed to take the continuum limit of Yang-Mills theory in the {x^,x'^} 
directions, and tune couplings to impose full Poincare invariance at finite transverse cut-off 
a. This procedure can be carried out using light-front quantisation about M = 0. Although 



this regime apparently cannot contain the Poincare-invariant theory at a = 0, numerical 
evidence for the existence of a renormalised trajectory was given, and has been extended to 
3 + 1 Yang-Mills ||13[ . In this paper, we present conclusive numerical evidence for the case 
of 2 + 1 Yang-Mills. 

For practical calculations, the remaining allowed operators in the action must be pared 
down to a finite number of independent parameters, and one must find some reasonable cri- 
teria to test Poincare invariance. We now develop these necessary approximations, following 
closely our previous work. 



A. Pure-glue sector 

To reduce the space of couplings to a finite dimension, we use various approximations: 

1. quadratic canonical momentum operator P"*", 

2. light-front momentum- independent couplings, 

3. transverse locality, and 

4. expansion in gauge-invariant powers of M. 



The reasoning behind them is described in more detail in Ref. ||13[. We only note here, 
that a poor choice of approximation will simply mean that we cannot get close to the 
renormalised trajectory, if it exists, and accuracy will suffer accordingly. The principle 
physical approximation. Item ^, is the 'colour-dielectric' expansion about M = 0, which is 
applied to the light-cone gauge-fixed Hamiltonian rather than the action. 

We have studied the light-cone Hamiltonian derived from the following SU{N) gauge- 
invariant action in the large- iV limit 

A = J dxUx' J2DaM{x')(D''M{x'))^ - V,. - ^Tr{F°^F,^} (2.1) 

where a G {0, 2} and 
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DaM{x^) = ((9„ + iA^{x^)) M(x^) - iM{x^)A^{x^ + a) , (2.2) 

is the tranvserse lattice covariant derivative. The 'potential' term is 

V^i = fi^ Tr[M{x^)M\x^)} + ^ Tt[m{x^)M\x^)M{x^)M\x^)} 

+ -^Tt[m{x^)M{x^ + a)M\x^ + a)M\x^)] + -^[Tr[M{x^)M\x^)]y . (2.3) 

In light-cone gauge A_ = and after eliminating A^ by its (constraint) equation of motion, 
the corresponding light-front Hamiltonian is 

P- = j:Jdx--^ tJj^x')^J^x')\ + ^ Tr J+(^^)^ Tr J^x') , +V^. (2.4) 

J+{x^) = i (^M(x^) d- M^(x^) + M^{x^ - a) H_ M{x^ - a)^ . (2.5) 

This is the most general Hamiltonian to order that obeys the other stated approx- 
imations. It can be light-front quantised and studied in a suitable Fock space at gen- 



eral momenta P+ and P^, as detailed in Ref. |]rT|. The eigenvalues of the exact Yang- 
Mills Hamiltonian yield the glueball masses M through the relativistic dispersion relation 
p- = (M2 + (pi)2)/2P+. 

B. Heavy Sources 

We introduce heavy sources (j){x^ , x~ , x^) on transverse sites. They are in the funda- 
mental representation and of mass p. We apply the same approximations that were made in 
the pure-glue sector, but here we expand to order all operators containing heavy-source 
fields, and work at leading non-trivial order in 1/p. The heavy-source action is Ag = A + A^ 
where 

A^ = j dx+rfa;-5](D„0)^D>-pV0-^Tr{P"^(a:i)P„/3(a:')W^(a;')} 

1i[M\x')F^^{x')M{x^)F^p{x' + a)] (2.6) 
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and 

W{x^) = (^M\x^)M{x^) + M{x^)M\x^)) . (2.7) 

Da = da + iAa{x^) is the usual covariant derivative for the plane x^}. After gauge fixing 
A_ = 0, eliminating in powers of M from its constraint equation, and discarding the 
higher orders in M, the Hamiltonian resulting from Ag which satisfies the approximations is 

^Trj^^W'j + ?^Tr{:^M(.')£^i^Mt(.>)j (2.8) 

with 

J+, = J+ + z</.H-0^ . (2.9) 

Like P^, Pg can be studied in a suitable Fock space. The eigenvalues of f^-P,", for co- 
moving heavy sources of velocity f^, are the usual excitation energies associated with the 
heavy-source potential 0]. If two sources are separated by na in the transverse direction x^ 
and by L in the longitudinal direction a;^, then a rotationally invariant string tension would 
imply that, for large separations, 

v^P~^aR, R = Va^n^ + (2.10) 

for the lowest eigenvalue. Demanding this rotational invariance, then comparing results 
at n = with L = 0, allows one to determine a in a dimensionful unit (we use G'^N) 
independent of a. This fixes the relative scale between x^ and x^, which will be needed 
for testing covariance. In practice it is relatively difficult to calculate the heavy source 
potential in the purely transverse direction. Consequently, we measure the string tension in 
this direction by compactifying space and calculating the winding mode spectrum. 

C. Poincare Invariance 

We test Poincare invariance of the theory by making measurements on eigenstates of 
P~ and P^ . It turns out that a rather simple set of tests suffices to obtain an accurate 

7 



estimate of the renormalised trajectory. One of the approximations made in arriving at P~ 
( p. 41 ) is transverse locahty. Therefore, it make sense to expand eigenvalues at fixed momenta 
(P+, P^) in powers of transverse momentum thus 

2P+P- = G^N [Ml + Ml a\P^f + Ml a\P^Y + ■ ■ ■) . (2.11) 

Note that G has dimensions of energy, and G'^N is held finite in the N —* oo limit.Q 
Mo, Ml, M2, ■ ■ ■ are dimensionless numbers which we calculate when diagonalising P~ . 
The simplest requirement of covariance is that 

Mla^G^N -1 = (2.12) 

This ensures isotropy of the speed of light. The dimensionless quantity a^G'^N has already 
been determined above from the scale setting procedure via the string tension. Further 
conditions come from higher order corrections in P^ in ( p.ll| ). In this work we will use only 
the condition ( |2.12| ) for the lowest-mass glueballs, together with conditions of rotational 
invariance in the heavy-source potential, to test the space of couplings of P^ . If our reasoning 
is correct and our approximations valid, we should find a well-defined trajectory on which 
the conditions ( |2.12| ) are accurately satisfied — in practice we introduce a test to quantify 



this. Moving along this trajectory should correspond to changing the spacing a. Eventually 
this would take us to the transverse continuum limit, but we will be prevented from reaching 
a = by the restriction /i^ > 0, a necessary condition for quantisation about M = 0. 

III. RESULTS 
A. charts 

It is convenient to form dimensionless versions of the other couplings 



^aG^ ^ 5^ as a — > 0, but since we do not approach a = we cannot use the continuum gauge 
coupling g. 
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2 \ 

The basic technique we follow is to search this space, with the test, for an approximation to 
a renormalised trajectory on which observables show enhancement of space-time symmetries 
violated by the cut-off a. The t^st is made up of variables to test isotropy of the speed 
of light in dispersion of low-lying glueballs, rotational invariance of the string tension, and 
rotational invariance of the potential at intermediate source separations. Since we can expect 
to do better with some variables than others, the weights are adjusted until we produce a 
sharp trajectory in coupling space where is minimized to roughly one per effective degree 
of freedom. In fact, altering the weights typically changes the sharpness of the trajectory 
and not its location. The optimum trajectory is tabulated in Table |. Full details on our 
computations are available at |17 



Figs. |I| and ^ show charts for a range of values of m vs. h and m vs. I2 near the 
renormalised trajectory.0 In each case the renormalised trajectory appears at the bottom 
of a well-defined and unique x^-valley, running from large to small m. The behaviour of 
the lattice spacing as one moves along the renormalised trajectory is shown in Fig. |[ As 
expected, the lattice spacing gradually decreases with but never becomes zero for > 0. 
The fluctuations are due mainly to the difficulty in establishing the scale ^/o^. Since the 
is stable and small over a range of lattice spacings, we will use the point of smallest lattice 
spacing (m = 0.044) to extract physical quantities. 

B. Rotational invariance 

The heavy-source potential is displayed in Fig. ^. It shows better restoration of spatial 
symmetry than previously obtained ||lTl. The potential in the continuum spatial direction 
x"^ is a fit to the lowest eigenvalue as a function of L of the form 



The behaviour of I3, which is always very large, is clarified in Ref. |11] 
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n 1 78 

V+P- = 0.15ALG^ N + 0.183 Vcm -—— . (3.2) 

Ij 

One must be careful when interpreting ( p.2| ) since the Coulomb potential in 2 + 1 dimensions 
is logarithmic. The form ( |3.2| ) should be appropriate except at the very smallest L, where 
Coulomb corrections are expected. The 1/L term is a universal correction expected on the 
grounds of models of flux-string oscillations ||18|. Universality implies that its coefficient 



should be invariant along the renormalised trajectory. In fact, we find that it varies slowly, 
a symptom that our approximation to the renormalised trajectory is not an exact scaling 
trajectory for this quantity and/or the form ( p. 21) is not sufficient to fit the potential. 

C. Glueballs 

The spectrum of glueballs in 2 + 1 dimensions can be classified by \J'\'^'^, where JT" 



is 5*0(2) spin, C is charge conjugation, and the parity V is spatial refiection 



Combinations of ±J' form parity doublets if states are Lorentz covariant. On the transverse 
lattice, there is enough symmetry to determine C, V and \J'\ mod 2. Additionally we can 
examine the shape of wavefunctions to help distinguish the spin of states. 

The lightest glueball is a O^"*"; its mass along the renormalised trajectory is shown in 
Fig. ^ The anisotropy of the speed of light in the 0+"'' dispersion is less than 3% for all the 
low points. For the point of smallest a, M.Q++ = 4.10(13)-\/o". Here, we have estimated 
a 2-3% error from extrapolations in DLCQ and Tamm-Dancoff cut-offs based on known 
analytic behaviour [^, and another 2-3% from systematic finite-a errors. The fractional 



finite-a error estimate is based upon deviations from the relativistic dispersion condition 
( p.l2| ). Figure ^ also shows the result of Teper, J^q++{N oo) = 4.065(55) who fit his 
finite-A^ data to a form A + B/N"^ in order to estimate the large-iV limit. Teper's large-A^ 
extrapolation and our independent direct calculation are in agreement .0 



^ We note that there are other finite-A^ lattice results which do not agree with Teper's. Recent 
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Fitting Teper's finite- data together with our large- result we find 

•^0- -4.118(13,^1^ + ^ (3.3) 



See Fig. |^. This gives a better estimate of the large-A^ limit than using one or the other 
data set alone. 

Although we have made improvements to the calculation of a, the dominant error in 
Fig. 1^ still comes from the fiuctuation of this quantity; in particular, the determination 
of o" from the longitudinal direction is a big source of error in determining the relative 
scales. This error becomes so severe for most heavier glueball states, which exhibit poor 
covariance, that an alternative method must be used for accurate results. To remove most 
of the error when dealing with heavier glueballs, we set Aio++/^/a to the large- value 
estimated in Eqn. ( p.3| ), then recalculated the renormalised trajectory with this constraint, 
id est we calculate mass ratios. To improve covariance in the lighter glueballs, at the expense 
of heavier states, we also restricted the to test only the dispersion of the lowest states in 
each charge-conjugation sector. The resulting mass ratios, at the point of lowest on the 
new renormalised trajectory, are shown in Table |l[ We also show the fit to the form 

^ (3.4, 



Mo++ m m 

including Teper's data. The convergence in is illustrated in Fig. |^. 



IV. CONCLUSIONS 

We have found that our improved transverse lattice calculations for 2 + 1 Yang-Mills 
in the large- limit are consistent with existing finite- A^ data from an independent lattice 
method. Although both make use of lattice regulators, the methods use different quantisa- 
tion procedures, elementary degrees of freedom, regulators, gauge fixing, and renormalisation 



Hamiltonian lattice calculations Q yield M.Q++ = 3.88(11)^0" for SU(3) compared to M.Q++ = 
4.329(41) in Ref. Q. 
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techniques. By combining the finite- and large- results, we have obtained accurate esti- 
mates of the lightest glueball masses and mass ratios for any in 2 + 1 dimensions. These 
should provide a useful benchmark for analytic studies. We are able to confirm that 0{1/N'^) 
corrections to the large- iV limit are typically small. The plots of mass versus could 
have had any shape, but they turn out to be almost straight and almost fiat. The same 
result also seems to be true in 3 + 1 dimensions [l^, though the data is less precise there. 



Recalling how the quark model explains the (OZI) suppression of corrections in most 
channels it would be interesting to know if constituent gluon models could provide an 
intuitive explanation of our finding. 

The lightest glueball mass in units of the string tension has recently been calculated 
analytically in 2 + 1 Yang-Mills for any Various attempts have also been made to 

obtain excited glueball mass ratios via extensions the ADS/CFT correspondence at large- 
A^ ||5|. These both use strong-coupling approximations and give qualitatively reasonable 
results. The transverse lattice method is also a strong-coupling one in a sense; it is a coarse 
lattice method. How can these methods give good results when the large- A^ approximation 
almost certainly introduces a phase transition in the strong-coupling/coarse-lattice regime 
The answer, we believe, lies in the fact that an exact renormalised trajectory, existing 
in an infinite-dimensional space of couplings, can avoid such transitions. Thus, it is only 
necessary to approximate the renormalised trajectory with a finite number of couplings to 
obtain results relevant to the continuum limit. One cannot extrapolate to weak coupling or 
fine lattice limits, but if one chooses the right variables and couplings, results can be obtained 
directly from the approximation to the renormalised trajectory far from these limits. 
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TABLES 



m 


h 


h 


h 




t2 




0.044 


-0.052 


-0.112 


680.2 


-0.661 


-0.691 


7.02 


0.089 


-0.087 


-0.109 


396.8 


-0.780 


-0.811 


7.85 


0.134 


-0.108 


-0.091 


3.221 


-0.896 


-0.876 


7.56 


0.180 


-0.147 


-0.107 


4.401 


-0.943 


-0.927 


6.55 


0.226 


-0.204 


-0.134 


178.5 


-1.098 


-1.167 


8.02 


0.2765 


-0.240 


-0.153 


5.48 


-0.989 


-1.138 


7.31 


0.3275 


-0.308 


-0.157 


6.01 


-1.181 


-1.340 


8.64 



TABLE I. The trajectory in coupling-constant space which minimises the t^st of covariance. 





M/Mo++ 


Fit coefficients 
C B 


A 


Q— 


1.35(5) 


-14.58(1.47) 2.983(191) 


1.349(6) 


2++ 


1.60(17) 


3.233(2.724) -1.144(856) 


1.743(51) 


o;- 


1.82(6) 


-5.839(7.488) 1.136(941) 


1.824(25) 


2-+ 


1.77(?) 


0.659(246) 


1.697(57) 


or 


1.28(?) 


0.770(399) 


1.520(38) 



TABLE IL Mass ratios for lightest glueball excited states, showing our N = oo measurement 
and fit coefficients including finite- data from Ref. [|| . The 2 ^ and 0^ states were not covari- 
ant-enough for reliable error estimates, and only Teper's finite- extrapolation is shown. 
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FIGURES 
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FIG. 1. Minimum for a given m and li. In the blank region to the right, tachyons appear 
in the spectrum. 




m 



FIG. 2. Minimum %^ for a given m and Z2. 
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FIG. 3. Variation of the transverse lattice spacing along the Lorentz trajectory. The fit is 
1.275m + 1.23. 
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FIG. 4. The heavy-source potential. Solid line is fit to potential for sources with x^-separation 



only; data points are values at one-link transverse separation and x^-separation Ly/G^N = 0, 2.5, 5. 
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FIG. 6. Fit to lowest glueball mass as a function of N. 
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